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We investigate competition between one- and two-dimensional topological excitations - phase slips 
and vortices - in formation of resistive states in quasi-two-dimensional superconductors in a wide 
temperature range below the mean-field transition temperature Tco- The widths w = 100 nm of our 
ultrathin NbN samples is substantially larger than the Ginzburg-Landau coherence length £ = 4nm 
and the fluctuation resistivity above Tco has a two-dimensional character. However, our data shows 
that the resistivity below Tco is produced by one-dimensional excitations, - thermally activated 
phase slip strips (PSSs) overlapping the sample cross-section. We also determine the scaling phase 
diagram, which shows that even in wider samples the PSS contribution dominates over vortices in a 
substantial region of current/temperature variations. Measuring the resistivity within seven orders 
of magnitude, we find that the quantum phase slips can only be essential below this level. 

PACS numbers: 74.78.Na, 74.40.+k 



The nature of the resistive state in superconductors at- 
tracts much attention from the physics community, since 
it involves fundamental phenomena and advanced con- 
cepts [l|, 0] such as the mechanisms of high-T c super- 
conductivity Q , thermal fluctuations [1, \Ei macroscopic 
quantum tunneling 0, 0, coherence S, topological 
excitations and phase disordering Resistive 

states are used in a number of quantum nanodevices, 
such as logic elements [l2[, ultra-sensitive detectors of 



radiation, single-photon counters, and nanocalorimeters 
3 14 1 . Understanding resistive states in nanoscale su- 
perconductors is critical for the advancement of funda- 
mental science and the development of novel applications. 

Phase slips and vortices are elementary topological ex- 
citations which create resistive states [l|, 0, [l(J. Wires 
with radius less than the coherence length £, or stripes 
with w < £ are one-dimensional (ID) superconductors. 
In ID-structures, the resistive state is produced by phase 
slips and is well described by the Langer, Ambegaokar, 
McCumber, and Halperin (LAMH) theory of thermally 
activated phase slips (TAPSs) 0,13,0,(1. At low enough 
temperatures, the quantum phase slips (QPSs) should be 
important but the magnitude of this effect and the 
characteristic resistance at the transition from TAPSs to 
QPSs are still under debate 

In 2D-superconductors the resistive state is formed 
by moving vortices. Above the Berezinskii-Kosterlitz- 
Thouless (BKT) transition temperature, Tc, there is a 
nonzero concentration of free vortices due to thermal un- 
binding of vortex-antivortex pairs (vap) Eon. Below 
Tc in an infinite 2D-superconductor, VAPs are tightly 
bound and only a significant bias current can unbind the 
pairs, resulting in a nonlinear flux flow resistance. In fi- 
nite size samples, free vortices can exist below Tc and 
produce a linear resistance at low bias currents [H| • It is 
commonly believed that the transition from ID phase slip 



excitations to 2D vortex physics takes place at w/£ ~ 1 
[ESQ EES- Howeve r, despite thorough studies of 
phase slip and vortex mechanisms, an investigation of 
their competition in quasi-two-dimensional superconduc- 
tors is long overdue. 

In this paper we study the resistive transition in NbN 
superconducting samples of thickness d = 4nm, width 
w = 100 nm, and a relatively short coherence length, 
£ ~ 4 nm, primarily due to the short electron mean free 
path. Thus, our samples are quasi-2D superconductors 
with a ratio w/£ = 25. We show that in these samples, 
ID phase slip strips (PSSs) are responsible for resistive 
state formation in a wide temperature range, extending 
far below the BKT transition. We also determine the 
phase diagram and show its scaling character. 

Our structures were fabricated from NbN supercon- 
ducting films, which were deposited on R-plane sapphire 
substrates by DC reactive magnetron sputterin g. D etails 
of the film deposition process are described in [18|. The 
samples were then patterned using direct electron beam 
lithography and reactive ion etching. Parameters for each 
of the samples are shown in Table I. 

Transport measurements were performed in a vacuum. 
4-point resistance measurements with RF-filtered leads 
were made with a lock-in amplifier for bias currents rang- 
ing from 10 to 300 nA. A 4-point DC measurement setup 
with RF-filtered leads was used for higher currents rang- 
ing from 500 nA to 1 /iA. In Fig.l we present the re- 
sistance vs. temperature obtained for both samples SI 
(Fig. l.a) and S2 (Fig. l.b) at a bias current of 10 nA. 

We begin by investigating the resistivity fluctuation 
region right above the mean field superconducting tran- 
sition Tco- The insets in Figs, l.a and l.b show the re- 
sistivity on a linear scale, just above the superconducting 
transition for both samples. Our analysis shows that the 
data obtained is well described by the Aslamazov-Larkin 
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Fitting the change in resistivity by 5p(T) = —<jalP%, we 
aim to determine Tco and pjv- The obtained parameters 
are presented in Table I. Note that we do not study the 
corrections very close to Tco and limit the range of data 
fitted to relatively small corrections, 8p/ pn < 0.1. Our 
fits in the insets of Fig. 1 show that in this temperature 
range the AL term in the form of Eq. 1 dominates over 
other possible contributions. 

Next, we investigate the resistivity below Tco m the 
framework of the LAMH theory for TAPSs. The LAMH 
expression for resistance at low bias currents is [l], 0, Q 
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where Rq — h/2e 2 ~ 13 kfi is the resistance quantum, 
£(T) = £(1 - T/T c) y 1/2 is the temperature dependent 
coherence length, and R sq = pN/d is the normal state 
resistance per square. 

In Fig. 1 we show fits of our R(T) data to LAMH 
formulas (Eqs. 2 and 3), which are in excellent agreement 
for both samples. In this fitting procedure we used only 
one fitting parameter £, which was found to be 4.0 run 
for both samples. This value is in a good agreement 
with £ « 0.6^hD/k B T ca = 3.4 nm, determined with the 
diffusion coefficient D ~0.5 cm 2 /s, which can be obtained 
from the resistivity of our films [201 ] . 

Measuring the linear resistance as a function of tem- 
perature in a range of seven orders of magnitude, we 
do not observe any deviation from LAMH theory, which 
is usually considered as a manifestation of QPSs. Al- 
though several groups reported observations of QPSs 
@> 0> 11] j the data and its interpretation remain con- 
troversial. While various phcnomcnological and micro- 
scopic models of QPS formation have been proposed, 
experimentalists appeal to the phcnomcnological theory 
suggested by Giordano [f|. In this approach, the QPS 
contribution to the resistivity can be described by Eq. 
2, where k B T is replaced by (8ak B /n)(Tco — T), with 
a a constant of the order of unity. Also, in this phe- 
nomenology an additional fitting parameter b was added 
into AF(T). According to Giordano, a is a universal 
constant, so the crossover from TAPS to QPSs should 
occur at a characteristic value of t cr = T /Tco indepen- 
dent on the sample cross-section and material parame- 
ters. The factor b is weakly dependent on material prop- 
erties and, therefore, the crossover should be observed 
at approximately the same value of r cr = R{t cr )/R^. 
In fact, for the first time the crossover was observed 
in In wires at r cr ~ 10~ 2 — 10 -3 [91; then in Sn wires 



at r cr ~ 10~ 3 - 10" 4 0, and recently in Al wires at 
r cr ~ 10~ 4 — 10~ 5 8]. All data shows a strong depen- 
dence of r cr on the wire diameter. Other recent mea- 
surements in MoGe wires do not show any trace of QPS 
for R/Rn spanning 11 orders of magnitude [20| . Our 
data supports the results of Ref. [20( and shows that 
QPSs could be important at lower temperatures and r cr 
values. 

We have also investigated nonlinear effects at bias 
currents ranging from 10 nA to 1 /iA. The current- 
voltage characteristics obtained at various temperatures 
are shown in Fig. 2. The solid lines represent the LAMH 
expression [lj], 



V{I,T) = 0.5 I T R LA MH(T)exp(I/I T ), 



(4) 



where I T = Aek B T/h ~ 0.013/zA/KT. To plot V(I,T) 
we used the resistance Rlamh (Eq. 2) investigated pre- 
viously at low currents. The dashed line in Fig. 2 repre- 
sents the characteristic current It- As shown in Fig. 2, 
the results of our measurements agree with the LAMH 
theory for currents up to ~ 1 /iA. Deviations at higher 
currents are most likely due to electron heating, which 
will be evaluated later. Note that, in thin-film quasi-2D 
superconductors, the heating can be localized into resis- 
tive domains, which are formed by the dissipative motion 
of vortices [9j; however, our data does not show such ef- 
fects of local heating. 

We have shown above that our data on resistivity mea- 
sured in ohmic and non-ohmic regimes are in excellent 
agreement with the LAMH theory, which has been de- 
veloped for quasi-lD superconductors. Why the LAMH 
theory turns out to be applicable to quasi-2D samples? 
Generally speaking, in superconductors with width wider 
than £, the order parameter can change along the width 
and, therefore, resistive domains of various geometries 
are possible. Eqs. 2 and 3 take into account only simple 
domains in the form of PSSs, i.e. strips with width equal 
to £ and a length equal to w (Fig. 3. a). Obviously, among 
other possible forms, the PSS has a minimal volume (see 
Fig. 3) and its generation requires a minimal energy, 
AF ~ (H 2 /87r)£W, where H c is the critical magnetic 
field. Because of the exponential dependence of Rlamh 
on AF, one can neglect domains other than PSSs. 

Overall, analysis of our data in the framework of the 
LAMH theory strongly suggests that the resistive state 
in our quasi-2D NbN samples is a result of ID PSS ex- 
citations. This observation is very interesting, because 
in such samples vortices are expected to dominate the 
resistive state. In what follows we support the statement 
above by modeling the competition between ID PSS ex- 
citations and vortices in quasi-2D superconductors. 

The vortex state critically depends on temperature 
with respect to the BKT transition temperature Tc, in- 
troduced for an infinitely wide film. Above Tc, there is a 
finite concentration of free vortices due to thermal fluctu- 
ations, while below Tc all vortices are tightly bound into 
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pairs. For disordered superconducting films, Tq is given 
by Tc/Tco = (1 + 0.54i? sg /i? )- 1 [21]. We estimate T c 
= 11.3 K for sample SI and T c = 11.6 K for sample S2. 
Since the obtained values of Tc are very close to Tco, 
and the temperature interval Tco — Tc is quite narrow 
compared to the entire resistive transition broadening in 
our samples, we limit our calculations to temperatures 
below Tc. In this temperature range, current-induced 
unbinding of VAPs can occur when the Lorentz force 
proportional to the bias current density exceeds the force 
of mutual attraction between the vortex and antivortex. 
The threshold current which unbinds the vortex and an- 



tivortex independent of film width is given by 23] 



I th ~2 K (T)I T , k(T) 



T c -T 
Tco - Tc 



(5) 



For I > I t h, the nonlinear resistance due to current- 
induced unbinding of vortex-antivortex pairs can be ap- 
proximated by 
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where Jo = [w / £(T)]It is a characteristic current propor- 
tional to the structure width. 

For I < I t h , the bias current is not sufficient to unbind 
VAPs, but a nonzero concentration of free vortices can 
still exist due to the finite size of the film [l6[ . These free 
vortices produce a linear resistance, 
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To investigate PSS and vortex contributions to the re- 
sistance, we use Eqs. 2, 4, 6, and 7. Taking into ac- 
count that the normal resistance may be represented as 
Rn — R S q{L/£)(£/w), we see that the relative contri- 
butions of PSS and vortices depend on four dimension- 
less parameters: I/I T , (T c - T)/(T C0 - Tc), to/f, and 
Rsq/Ro- In Fig. 4 we present a phase diagram of the 
resistive states in I /It and [Tc — T)/(Tco — Tc) coor- 
dinates for samples with w/£, — 25 (Fig. 4. a) and w/£ 
— 75 (Fig. 4.b). The value of R sq /Ro was taken to be 
0.05, which corresponds to our NbN samples. The solid 
lines in Fig. 4 represent the boundary between PSS and 
vortex mechanisms, described by Eqs. 2 and 4 and Eqs. 
6 and 7 correspondingly. The solid gray lines represent 
It which separate the linear and nonlinear mechanisms 
related to PSSs. The nonlinear vortex mechanism domi- 
nates over others at currents, which exceed the threshold 
current given by Eq. 5. The hashed area represents a re- 
gion of strong electron heating, where the change in elec- 
tron temperature 59 ~ j 2 pT e _ p h/C e rises to (Tc — T). 
In these evaluations we have used values for the elec- 
tron heat capacity, C e , and the electron-phonon relax- 
ation time T e -ph from Ref. 23j. Finally, the dashed 
segments of the boundaries in both phase diagrams rep- 
resent a narrow region very close enough to Tco, where 



the interaction between phase slips and nonequilibrium 
phenomena become important and, therefore, the LAMH 
theory is not applicable [l|. 

As the sample width increases, the relative contribu- 
tion of PSSs to the resistivity decreases. The phase dia- 
gram for samples with w/£ — 75 in Fig. 4.b shows that, 
at currents below I t h, the PSS mechanism competes with 
free vortices generated below Tc due to finite size effects. 
While the contributions from both mechanisms decrease 
with increasing sample width, the PSS formation is more 
sensitive to the parameter w/£ because of an exponen- 
tial dependence of the PSS resistance (Eq. 2) on AF(T), 
which, in turn, is proportional to w/£. Fig. 5 shows that 
the characteristic crossover temperature T* separating 
the linear PSS and vortex mechanisms strongly depends 
on the sample width and changes very weakly with R sq . 
Therefore, the obtained diagram has a scaling character 
in the dimensionless coordinates used. 

Reviewing other experiments with NbN, let us note 
that the complex investigations of resistivity in 1 (im 
films [24j . demonstrate only effects of vortices. This re- 
sult is in agreement with our modeling, which also does 
not show any measurable contribution of PSSs in the 
samples with — 250. In [25| samples with geometry 
similar to our device S2 and Tco=4.8 K have been inves- 
tigated at rather high currents, ~ 1 /iA. The data of Ref. 
25 have been found to be well described by the vortex- 
induced resistivity. This observation is also in qualitative 
agreement with our phase diagram, while essential effects 
of electron heating may also be expected. According to 
our modeling, samples with w/£ in the range of 50 - 100 
(0.2 -0.4 /j,m for NbN) would be preferable for further ex- 
perimental investigations of the boundaries between ID 
and 2D excitations. 

In conclusion, we have investigated the competition be- 
tween the PSS and vortex excitations in the formation of 
resistive states in quasi-2D superconductors. Our data is 
in excellent agreement with the LAMH theory of TAPSs 
and does not show any evidence of QPS contribution. 
The data and modeling show that the ID PSS mecha- 
nism dominates over the substantial range of parameters 
(see Fig. 4). Essentially, the creation of a PSS over- 
lapping the sample cross-section requires a significantly 
higher energy than the creation of a vortex or VAP. How- 
ever, the PSS mechanism turns out to be more effective 
in producing resistance. In clean superconductors with 
large values of £, the PSS excitations will prevail over 
vortices in wide samples, because the phase diagram de- 
pends only on the parameter 
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TABLE I: Sample parameters: the sample thickness (d), length (L), and width (w); the mean field transition temperature 
(Tco), the normal resistivity (pjv), and the Ginzburg-Landau coherence length (£). Tco and pjv were determined from the AL 
conductivity above Tco , £ is found from the LAMH fit below Too- 



Sample 


d 


L 


w 


Tco 


Pjv 






nm 


pm 


nm 


K 




nm 


SI 


4 ±0.2 


500 


100 ±5 


11.6 


264 


4.0 


S2 


4 ±0.2 


30 


100 ±5 


11.9 


252 


4.0 
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FIG. 1: Normalized resistance R/Rn vs. T for the samples 
SI (a) and S2 (b), the solid lines represent fits by the LAMH 
theory. Insets show the resistance above Too, and solid lines 
are the AL fits. 



FIG. 4: The phase diagrams for samples with w/£ of 25 and 
75. Hashed pattern represents the area of significant electron 
heating. 



FIG. 2: Current-voltage characteristics of the sample SI at T 
= 9.7 K - 10.6 K. The solid lines are fits by the LAMH theory 
(Eq. 4). The dashed line shows the characteristic current It- 



FIG. 5: The crossover temperature T* as a function of sample 
width for samples with different sheet resistances. 



FIG. 3: Among other possible forms of resistive domains 
in quasi-2D superconductors, the PSS stripe (a) with width 
equal to £ and a length equal to w has a minimal volume. 
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